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Abstract 

There is evidence for a 130 GeV 7-ray line at the Galactic Center in the Fermi Large 
Area Telescope data. Dark matter candidates that explain this feature should also annihilate 
to Standard Model particles, resulting in a continuous spectrum of photons. To study this 
continuum, we analyze the Fermi data down to 5 GeV, restricted to the inner 3° of the 
Galaxy. We place a strong bound on the ratio of continuum photons to monochromatic line 
photons that is independent of uncertainties in the dark matter density profile. The derived 
constraints exclude neutralino dark matter as an explanation for the line. 
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I. 



INTRODUCTION 



Astrophysical searches for dark matter are a critical component of the experimental effort to 
explore the dark sector ([U 12]). The strategy is to look for observable products of dark matter 
annihilation, both in the Milky Way halo and beyond. A variety of experiments are currently 
underway (see [3] for a review) , exploring many different annihilation products over a broad range 
of energies and target spatial regions. A particularly important potential astrophysical signal is a 
monochromatic 7-ray line, with energy corresponding to the mass of the annihilating dark matter. 
The discovery of such a line would be a "smoking gun" signature of dark matter because no 
astrophysical backgrounds are known to generate a peak in the 7 spectrum. 

The Fermi Large Area Telescope (LAT) is currently taking data, and continually improving its 
sensitivity to features in the 7 spectrum that could be signatures of dark matter [IHS] . A recent 
analysis of the Fermi 7 spectrum from 20-200 GeV has found preliminary evidence for a sharp 
feature around the Galactic Center corresponding to E-y ~ 130 GeV |10l lllj . The analysis by |10j 
finds 3.3 a evidence (including the look-elsewhere effect) for such a line in a region of the sky 
that extends roughly 15° above and below the Galactic plane. Based on an observation of ~ 50 
photons, the tentative signal corresponds to dark matter with best fit mass m x = 129.8 GeV and 
annihilation cross section a~~v ~ 1.27 x 10~ 27 cm 3 /s, if the dark matter density follows an Einasto 
profile. 

Following [10], the authors of |12t IT3] have confirmed the presence of the 130 GeV line and 
have strengthened the case that this excess could be due to dark matter annihilations. The au- 
thors showed that the signal is concentrated in a ~ 3° radius region about the Galactic Center. 
Additionally, they obtain better fits if the chosen region is off-set from the center of the galaxy by 
0(1°), and if the signal consists of a pair of lines at 111 and 129 GeV |13j . 

The origin of the 130 GeV feature in the Fermi spectrum remains a mystery. At present, the 
official search by the Fermi collaboration has found no evidence for monochromatic 7 lines [3] . For 
a 130 GeV dark matter mass, this analysis sets an upper limit of l.OxlO - 27 cm 3 /s for an Einasto 
profile, which is in tension with the purported signal. However, the two analyses are fundamentally 
different in their approaches. The Fermi collaboration searched for lines using the all-sky 7-ray 
maps; the analysis was done for |6| > 10° plus a 20° x 20° square at the Galactic Center using the 
Pass 6 data. In contrast, the search regions in |10j were defined to optimize the significance of a 
dark matter signal and the analysis was performed using the Pass 7 data. 

Our goal here is to explore constraints on the properties of dark matter that can be obtained 
from the Fermi data under the assumption that the 130 GeV line is due to annihilating dark 
matter. In a large class of weakly interacting dark matter models, one expects the annihilation 
of dark matter into a pair of photons to arise from loops of states that carry electroweak charges. 
One consequence is that there is a non-zero annihilation rate into 7Z and/or 7/1. If the particles 
in the loop are lighter than the dark matter, annihilation into these states can dominate and a 
continuum spectrum of photons results from the subsequent decay of the annihilation products. 
This is the case for the neutralino of the minimal supersymmetric Standard Model (MSSM) [2], as 
well as non-supersymmetric extensions of the Standard Model, such as Universal Extra Dimensions 
|14[ I15j . Little Higgs models |16(. I17j . and more generic weakly interacting (WIMP) dark matter 
models [EM]. 

In this work, we constrain the ratio of the number of continuum photons to the number of 
photons responsible for explaining the 130 GeV 7 line. The expression for the photon flux is given 



by 




(1) 
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^■rnax 


= 100 


DATA.QUAL 


= 1 


LAT.CONFIG 


= 1 


ABS(ROCK_ANGLE) 


< 52 



TABLE I: Criteria used to filter events. 



where <r ann is the total dark matter annihilation cross section into final states /, dn^/dE is the 
differential photon spectrum resulting from each /, Brj is the branching fraction into /, £(9) is 
the line of sight (LOS) as measured at an angle 9 from the Galactic Center, and p(£) is the dark 
matter density along the LOS. Constraints on the ratio of photons from different final states are 
particularly powerful because dependencies on the dark matter density profile cancel out. As a 
result, the constraints derived in this paper are independent of astrophysical uncertainties. 

We present two versions of the constraint on the dark matter continuum contribution. The 
first requires that the continuum photons not supersaturate the observed Fermi data. The second 
involves a log likelihood fit of the shape of the dark matter and background spectra to the Fermi 
data. The second constraint is much stronger, but relies on the assumption that the astrophysical 
background is well-described by a single power-law. Our derived constraints exclude the neutralino 
explanation for the 130 GeV 7 line. 

We begin in Sec. [II] by discussing the Fermi data used in this analysis and the details of our 
fitting procedure. In Sec. |III[ we provide the constraints on the ratio of continuum photons from 



annihilations to W + W and Z°Z° to the number of photons that yield the 7 line. In Sec. IV 



these constraints are applied to neutralino dark matter. We summarize our conclusions in Sec. [V} 
There are three appendices to the article that provide the counts per bin extracted from the public 
Fermi data used in this analysis (Appendix [A]), discuss our computation of the Fermi Instrument 
Response Function for energy dispersion (Appendix [B]) , and provide constraints for dark matter 
annihilations to the additional final states bb, t + t~, and (Appendix [c| . 



II. A 130 GEV GAMMA LINE 



In this section, we describe the Fermi data used in our analysis and the methodology employed 
in our log likelihood fitting procedure. As an application of our statistical procedure, we consider 
the best fit for a pair of lines, which is expected if the dark matter annihilates to both 77 and 7Z . 



A. The Fermi Data 



The Fermi LAT 3.7 year public data that is used in this analysis includes the direction of ar- 
rival and the reconstructed energy for each measured photon in the Pass 7_Version 6 release. 1 
A standard zenith angle cut is applied to exclude events at angles greater than 100°. The live- 
time/exposure is computed as recommended for diffuse analyses: 2 specifically, we use option 2 with 
the criteria specified in Tab. |Ij The DATA_QUAL filter excludes time periods when the quality 
of the data is compromised, e.g. solar flares, the LAT.CONFIG flag includes data from when the 
LAT is in nominal science configuration, and the ABS(ROCK_ANGLE) cut excludes periods when 



1 http://heasarc.gsfc.nasa.gov/FTP/fermi/data/lat/weekly/p7v6 

2 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/Cicerone_Likelihood/Exposure.html 
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the Earth limb is in the field of view. We use the class of events designated ULTRACLEAN, which 
have a lower effective area but also a lower background than the SOURCE class. 

Following |13| . we restrict our analysis to the inner 3° radius region around the Galactic Center 
and neglect possible enhancements from an offset along the plane [13] . Unless explicitly stated, 
all results use data where the area within 1 degree of the Galactic Center is masked to reduce 
background contributions. We restrict to the energy range 5-200 GeV to minimize uncertainties due 
to the point spread function (PSF). The Fermi LAT is designed to measure photons from around 
20 MeV to many hundred GeV. The PSF, which encodes the uncertainty in the reconstructed 
position in the sky, starts to grow rapidly below a GeV. Specifically, the 68% containment radius 
of the PSF is about 0.9° at 1 GeV and decreases with energy, approaching ~ 0.2° at high energies. 

Appendix |A| provides the counts per bin for the relevant region of the sky when the inner degree 
is both masked and unmasked. The photon counts are given for iVbins = 128 energy bins from 
5.1-198 GeV. 



B. Fitting The Data 

For concreteness, we assume that the signal arises from a WIMP of mass m x annihilating into 
77 and/or 7Z , thereby producing at most two lines in the photon spectrum at energies 

The WIMP may also annihilate into final states {e.g., W + W~ , Z° Z°, bb, t + t~ , n + fi~, etc.) whose 
decay products shower and hadronize to produce a continuum photon contribution. Assuming that 
the background is a falling power-law parametrized by a, (3, the observed photon spectrum expected 
from this model is 



cP(E) = C E a(E) 




+ N^D(E, £ 77 ) + N JZ D(E, E 7Z ) + (E, m 



(3) 



where iV 77 , iV 7z , and iV ann are the normalizations of the separate signal components. The function 
D(E, E true ) is the energy dispersion about the true signal energy and is derived using the Fermi 
Instrument Response Function (IRF) obtained from the publicly available Science Tools 3 — see 
Appendix [B] for a detailed discussion. The normalized differential distributions for different anni- 
hilation final states, denoted dn 7 /d-E, are obtained using Pythia version 8.165 [22] to generate the 
spectra. Cea is a corrective factor that accounts for the change in effective area in the 3° region 
about the Galactic Center, as a function of energy. 

For Poisson-distributed data, the best fit values of the parameters a, (3, iV 77 , -/V 7z , and iV ann are 
obtained by maximizing the likelihood function 

N hins 

]n£(a,P,Ny 1 ,Nry z ,N axm )= y~] n& • hi<ft fc - <fa - lnn fc !, (4) 

k=i 

where is the observed photon count and <pk = Jc;™ ax 4>(E)dE for the fc th bin spanning 

min 

[-^min> ^max] • The confidence region about the maximum likelihood, ln£ max , is determined by 

ln£ > ln£ max — Aln£, (5) 



3 http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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where 2 A ln£ = A% 2 and the number of degrees of freedom (d.o.f.) is the number of fit parameters. 
In Fig. [2} the 1, 2, and 3 a contours refer to Aln£ = 1.76,4.01, and 7.08 (3 d.o.f.). In Fig. [5} the 
1, 2, and 3 a contours refer to Aln£ = 2.36,4.86, and 8.13 (4 d.o.f.). 

Next, we use this statistical procedure to show that the photon spectrum in the region of interest 
is consistent with the presence of a photon line. For now, we assume that the photon continuum 
does not contribute to the signal, reserving the case where iV ann > for the next section. Scanning 
over m x and 

# 7z / 77 = arctan (6) 

while maximizing over a, 0, and iV 77 , we find that the best fit point corresponds to 

{m x /GeV,a,/3,N 77 ,# 7z/77 } max = {130,2.67,0.88,30.3,0} (unmasked); 

(7) 

{m x /GeV,a,/3, N 77 ,# 7z/77 } max = {130,2.62,0.80,31.6,0} (masked), 

where masked (unmasked) refers to removing (including) data within 1 degree of the Galactic 
Center. The significance of the best fit point relative to the null model (power-law background) is 
5.5 a for both the masked and unmasked cases, not including look-elsewhere. 4 Masking a 1 degree 
radius circle around the Galactic Center has little effect on the best fit dark matter parameters, 
though it prefers more shallow power-law backgrounds. From this point onwards, we will only 
consider the masked data. 

Figure [T] shows the spectrum of photon counts in the region of interest. The solid red line 
corresponds to the best fit model in Eq. (Frj) obtained by maximizing the likelihood function over 



T 




E y [GeV] 

FIG. 1: Photon counts within 3° degrees of the Galactic Center with the inner degree masked. The solid 
red line shows the best fit model given in Eq. 0, assuming no continuum contribution. The dashed 
black line shows the continuum spectrum for a 130 GeV dark matter annihilating into W + W~ (arbitrary 
normalization); the spectrum for Z°Z° is indistinguishable. 



4 The best fit null model is {a,/3} n uii = {2.65,0.95} for the unmasked case and {a,/3} n uii = {2.58,0.87} for the 
masked case. 
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FIG. 2: Regions of 1, 2, and 3 a significance (filled contours) for # 7z / 77 = arctan N lz /N 11 as a function of 
mass for the case N ann = 0. The solid lines are contours of JV 77 + N lz . The best fit point, marked with a 
white cross at m x — 130 GeV and # 7z / 77 = 0, is given in Eq. (|7j). This figure was made using the masked 
data. The analogous plot for the unmasked case is qualitatively the same. 

the energy range from 5-200 GeV. The spectrum is well-characterized by a single falling power-law 
and a peak at 130 GeV comprised of ~30 photons. 

Figure [2] shows the 1, 2, and 3 a contours for points in the # 7z / 77 — m x plane. The best fit 
point is marked by the white "X." There is a clear symmetry in the significance contours, with 
regions about 130 and 145 GeV each within 1 a of the best fit model. For the case of a 145 GeV 
dark matter, all the photons in the 130 GeV line are due to dark matter annihilation to 7Z . 

Both [13] and [23] note that the presence of two lines at ~ 130 and 115 GeV is a slightly better 
fit to the data; we reproduce the results in |23] when redoing our analysis with the data from jlOj . 
However, for the region of interest analyzed in this work, the data prefer a single line corresponding 
to a DM mass of either 130 or 145 GeV, although two lines are consistent within 1 a. 

III. CONTINUUM CONSTRAINT 

Now that we have demonstrated the presence of a line (or pair of lines) in the Fermi data 
at the Galactic Center, we explore correlated photon signals that are important for a wide class 
of models. If the dark matter annihilates to Standard Model particles beyond 77 and/or 7Z , 
additional photons are produced in the decay of these states. The Fermi data has been used 
to place constraints on the resulting inclusive photon spectrum using a variety of methods [U- 
El [2"4"1 [2"5~] . Currently, the strongest of these bounds comes from a search for 7 rays originating from 
10 dwarf galaxies in the Milky Way over a 24 month period [6]. For annihilations to W bosons 
and m x ~ 130 GeV the 95% confidence bound is 

c ww v < 10~ 25 cm 3 /s (Fermi dwarf Galaxy constraint). (8) 

The constraint from stacked dwarfs rules out large regions of parameter space for models that can 
potentially explain the 130 GeV 7 line. In the MSSM, for example, assuming the neutralino makes 
up all the dark matter, this constraint rules out all models with \a lz v + cr 77 f > O (l0~ 28 cm 3 /s), 
except for wino-bino mixed neutralinos. 
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Here, we derive a constraint on the ratio of the number of continuum photons to the number of 
photons in the peak using the data from the Galactic Center. Specifically, we constrain the ratio 



R th = ^ann ^ (g) 

2 (Tsyy ~\~ C^'Yt' 



derived for a given theory, 5 by comparing it with the associated quantity obtained from observation 

pob _ 1 ^ann ,-. n s 

= n 7 N +N ' ( } 

where o" ann is the total dark matter annihilation cross section, and cr 77 (c 7z ) is the annihilation 
cross section to 77 (7Z ). -/V ann is the number of photons in the continuum spectrum that results 
from the process that dominates cr ann , and iV 77 (N^ z ) are the number of photons in the peak(s) 
attributed to dark matter annihilations to 77 (7Z ). n2nn is the total number of photons per 
annihilation in the considered energy range. We do not include this factor in the definition of i? th 
because it depends on the energy range of interest and will be different for the "supersaturation" 



constraint (Sec. Ill A) and the "shape" constraint (Sec. IIIB), discussed below. 

In the following subsection, we constrain R oh by requiring that the continuum photons do 
not supersaturate the data. The result is model-independent, and only depends on the final 
state annihilation products. Because the bound is a ratio, it is independent of any astrophysical 
uncertainties and applies to scenarios where the annihilating particle is a subdominant component 
of the dark matter [SB]. We then go on to present an even stronger bound using the full shape 
information of the continuum spectra. Although this constraint is significantly stronger, it does 
depend on whether a single power-law explains the background from 5-200 GeV. 

In what follows, we focus on dark matter annihilations to W + W~ and Z°Z°. Shape constraints 
for final states bb, t + t~, and /z + 'fjT are given in Appendix [cj 



A. Constraint from Supersaturation 



In this subsection, we derive a constraint on R ob as defined in Eq. (10) that is independent 
of any background model assumptions. This constraint arises from the fact that the continuum 
contribution should not supersaturate the data. It is conservative in that it assumes that the 
entirety of the photon spectrum is due to signal, with no background contribution. 

To obtain the optimal supersaturation exclusion, we must select an energy bin where the number 
of continuum photons peaks relative to the background. The spectral index a of the 7 continuum 
background is expected to follow that of the proton spectrum a p from 5-100 GeV (see [27J and 
references therein). For concreteness, we use the measurement from the PAMELA collaboration 
in the range E p = 30 - 80 GeV: a p = 2.801 ± 0.007 (stat) ± 0.002 (syst) [28]. Note that this is a 
local measurement and may differ at the Galactic Center. It does, however, provide an independent 
determination of the spectral index and gives a reference point from which to select an optimal 
energy bin. 

In Fig. |3| we show multiplied by the differential continuum photon spectrum as a function 
of E-y, where E^ is the photon energy, for W + W~ , Z° Z° [black, solid], bb [green, dashed], t + t~ 
[blue, dotted], and [red, dot-dashed] final states. The location of this peak is approximately 



5 When R th is 0(1), one must be careful to include the contribution to the continuum spectrum from \X 1%° 1 
which has a slightly different shape than \X ~^ Z a Z° from kinematic effects. We neglect this subtlety for the 
constraints presented in this paper. 
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FIG. 3: Continuum photon spectra x E^ 8 as a function of E n for dark matter annihilating into 
W + W~ , Z° Z° [black, solid], bb [green, dashed], t + t~ [blue, dotted], and [i + ^ [red, dot-dashed] . The 
dark matter mass is taken to be 130 GeV. 



at 15 GeV for W + W~ , indicating where the continuum spectrum should peak over the power-law 
background for this final state. 

The shapes of Z° Z°, and bb fall off sufficiently fast as E 1 —t m x . However, for 

t + t~ and the continuum spectrum can make a non-trivial contribution at ~ m x . For 

the supersaturation constraint, the number of photons in the peak must be independent of the 
continuum, and it is therefore important that the continuum contribution to energies near m x be 
negligible. This is the case for W + W~ , Z° Z°, and bb, but not for t + t~ and fi + fi~. Therefore, we 
will not compute the supersaturation constraint for the leptonic final states. 6 Due to the similarity 
in shape between W + W~ and bb, we only present the supersaturation constraint for the former. 

In Fig. [4j we show the 95% C.L. exclusion region for the supersaturation analysis. The "optimal 
bin" is 10 — 20 GeV and contains 1104 photons. To determine the exclusion region, we assume 
that the number of events in this bin are Poisson distributed and marginalize over the number of 
photons attributed to the 7-line signal (using a fit to the peak for data from 80-200 GeV). For 
reference, there are about 29 photons contributing to the 7-line at m x = 130 GeV and 145 GeV 
- the results are qualitatively similar to those shown in Fig. [2] To account for the change in 
effective area between the two bins, we apply a multiplicative correction of 1.06 to the number of 
peak photons; this is computed by taking the ratio of the effective area at 130 GeV and 15 GeV. 
Pythia is used to determine the average number of photons per annihilation, n2nn, for the process 
XX -> W+W- in the range 10-20 GeV. For m x = 130 GeV (145 GeV), nj nn = 0.67 (0.78). 

Requiring that dark matter annihilations to W + W~ do not supersaturate the data at the 95% 
C.L. constrains R ob > 94 (77) for m x = 130 GeV (145 GeV). Note that the constraint weakens at 
the edges of the considered mass range; at these masses, a line does not provide a good fit to the 
data and therefore the best fit value of iV 77 + N~ jz goes to zero, causing i? ob to become larger. To 
good approximation, the limit on R oh is the same for the bb final state. 



Sec. IV will show that there is a lower bound of R th > 200 in the MSSM for annihilations to 



The fact that r + r and fi + n final states contribute photons at energies near the dark matter mass also affects the 
shape constraints and accounts for the differences in the excl usion curves for W + W, Z°Z°, and bb (see Figs. [5] 



and 9| and r r and [jl [A final states (see Figs. fO and 11\ 
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FIG. 4: The 95% C.L. excluded region for R oh , as defined in Eq. (flOb, versus m x assuming annihilation 
into W + W~ , Z° Z° for the supersaturation analyses using the masked data set. The plotted mass range 
corresponds to the 2 a best fit region. For comparison, R^ ino — 200 and -Rniggsmo — ^00. Pure wino and 
Higgsino dark matter are clearly excluded, as discussed in Sec. |IV| 



W + W and Z°Z°. Clearly, the supersaturation constraint robustly rules out this entire parameter 
space. 

B. Constraint Utilizing Shape Information 



In this section, we present a complementary bound on R oh that utilizes the shape of the con- 
tinuum spectrum. The ratio iV ann / (-/V 77 + N lz ) is constrained by performing a log likelihood fit as 

and 



described in Sec. II B For a given value of iV a 



/{N^ + N lz ) 



"X: 



we marginalize over a, j3, 



N^, and N~ z . This analysis is more constraining than the supersaturation results of Sec. Ill A 



but depends on the assumption that the 7 ray background is described by a single power law from 
5-200 GeV. 

The best fit point is the same as in Eq. ([7]), with -/V ann = 0. The fact that the fit prefers 
no annihilation to W + W~ is not surprising. Figure [l] shows that a single power law provides a 
remarkably good fit to the data between 5-100 GeV. The filled contours in the left panel of Fig. [5] 
show the 1, 2, and 3 a confidence regions about the best fit point. The black solid lines denote 
contours of iV 77 + iV 7z . There is some room for a non-zero annihilation contribution. For these 
cases, the continuum spectrum explains the data below ~15-20 GeV and the power law background 
becomes important at larger energies. Typically, the best fit power law is shallower when -/V ann > 
than when -/V ann = 0. 

The 2 a confidence region for iV aiin /(jV 77 + N lz ) can be converted into a bound on R oh by 
multiplying by l/n^nn integrated over the appropriate energy range. The result is given on the 
right in Fig. [HJ which shows the region excluded at 95% C.L. for R oh . The maximum allowed value 
is i?max — 10 f° r a mass of 129 GeV. The entire range of R oh is excluded outside the plotted range 
for m x because these masses do not provide a good fit to the data. 

Electrons and positrons produced by dark matter annihilation can give additional contributions 
to the continuum from inverse Compton scattering (ICS) of the interstellar radiation field [29J. Ne- 
glecting this contribution is conservative for the supersaturation constraint, but one might wonder 
if the addition of ICS photons could improve the spectral correspondence between the model and 
the data for the shape constraint, hence weakening the limits. 
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The details of the ICS contribution depend on the diffusion and transport of charged particles 
in the inner Galaxy; however, for the W + W~ , Z°Z°, and bb annihilation channels, the ICS contri- 
bution is always subdominant. To demonstrate this point, consider the case where all the energy 
of the electrons and positrons is converted into upscattered photons via ICS (as opposed to other 
processes such as synchrotron), and furthermore where spatial diffusion can be neglected, with all 
scatterings occurring at the point of annihilation. These approximations likely overestimate the 
true ICS contribution; any significant energy density in magnetic fields will reduce the power in 
ICS photons in favor of synchrotron, and since we are examining the region where the line sig- 
nal is concentrated, diffusion of electrons and positrons after annihilation should only dilute the 
corresponding ICS signal. 

We obtain the electron and positron spectra produced by DM annihilation using Pythia as 
described previously, and then compute the spectrum of photons produced by repeated ICS (equiv- 
alently, the ICS photons produced by scattering on the steady-state electron+positron spectrum), 
under these simplifying approximations. Scattering of ~ 100 GeV photons on starlight approaches 
the Klein-Nishina regime, where the energy of the upscattered photons is comparable to the initial 
electron energy; consequently we use the full Klein-Nishina cross section and take into account 
that in this regime the electrons may lose a significant fraction of their energy in a single scattering 
(see [30j for formulae and discussion). We employ the model for the interstellar radiation field at 
the Galactic Center used in GALPROP version 50p [31j ; while there could be other radiation fields 
present, in order to dominate the ICS losses they would need to have an energy density exceeding 
that of the CMB or the starlight photons in the inner Galaxy, which seems implausible. 

As shown in Fig. [6j for dark matter annihilations to W + W~ , the ICS spectrum is everywhere 
subdominant and adding it makes no significant difference to the continuum spectrum. The same 
conclusion holds true for Z°Z° and bb. As a result, it is not necessary to include the ICS contri- 
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FIG. 5: (Left) The 1, 2, and 3 a confidence regions (filled contours) for N ann /(N 71 + N lz ) as a function of 
mass for dark matter annihilation to W + W~ . The ratio N lz /N 11 is allowed to freely vary for each point in 
the grid. The solid black lines are the contours for N 11 + N lz . The best fit point is marked with a cross at 
m x — 130 GeV, # 7z / 77 = 0, and N ann = 0. The confidence regions for Z°Z° and bb are similar. (Right) The 
shape analysis constraint. The shaded region corresponds to parameters where the fit is 2 <r or worse with 
respect to the best fit point. This constraint is O(10) stronger than the supersaturation constraint shown 
in Fig. H The constraints for bb, t + t~, and /i + ^t~ final states are provided in Appendix [cj 
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FIG. 6: We show the spectrum of photons produced by the decays of W + W~ from a single DM annihilation 
[orange, dotted], an estimate of the spectrum resulting from inverse Compton scattering of the associated 
electrons and positrons on the microwave, infrared and optical photons in the inner Galaxy [gold, dashed], 
and their sum [black, solid]. The dark matter mass is 130 GeV. We see that for this final state, the photons 
from ICS are a negligible component of the full continuum photon spectrum in the energy range relevant to 
this study. 



bution in our fitting procedure. 

At low energies, the ICS contribution for leptonic final states computed using this simple esti- 
mate dominates over the directly produced photons. We have explicitly checked the effect on our 
likelihood analysis. The inclusion of the (likely overestimated) ICS contribution described above 
only tightens the limits; it makes the constraints more stringent while not shifting the best-fit 
points. Thus, the bounds we present on leptonic final states are conservative. Because the best-fit 
value for iV ann is non-zero for these final states, a more careful treatment of the ICS contribution 
may be worth future study. 



IV. IMPLICATIONS FOR NEUTRALINO DARK MATTER 

Now that we have explored the constraints on the dark matter scenario with both line and 
continuum contributions, this section will be devoted to the implications for MSSM neutralino 
dark matter. The MSSM neutralino is a linear combination of bino B, wino W, up-type Higgsino 
H u , and down- type Higgsino Hd- Its properties are controlled by four parameters that determine 
the mass matrix: the bino mass M\, the wino mass M2, the Higgsino mass ti, and tan/3. The 
DarkSUSY program version 5.0.5 |32] is used to compute all relevant annihilation cross sections 
|33H36j . The conventions for the neutralino mass matrix are the same as in |37| : we do not impose 
any priors on the relic density. 

To understand the range of possible values for cj 77 , a lz , and cr a nn in the MSSM, we begin by 
analyzing the pure gaugino eigenstates for m x = 130 GeV. In the decoupling limit (tha —> 00) with 
heavy sfermions, the pure bino is inert and only pure wino and Higgsino states have non-negligible 
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annihilation cross sections: 

a^v ~ 2.5 x 1CT 27 cm 3 /s 
Wino: a lz v ~ 1.4 x 1(T 26 cm 3 /s => # th = 210; (11) 

O-annV ^ (J w t> ~ 4.0 X 10" 24 Cm 3 /s 

and 

(J 77 v ~ 1.1 x 10~ 28 cm 3 /s 
Higgsino: <j 1z v ~ 3.7 x 10~ 28 cm 3 /s => i? th = 710. (12) 

O-annV - + ^ 4.2 X 10~ 25 Cm 3 /s 



For these two cases, the cross sections to 77 and 7Z are dominated by loops involving the 
charginos. Annihilations to W + W~ are due to i-channel chargino exchange. The Higgsino also has 
a non-trivial cross section to Z°Z° due to the presence of a light Higgsino-like second neutralino. 
The pure Higgsino and pure wino are clearly ruled out by the constraints presented in Fig. [4} 

To explore the case where the neutralino is a non-trivial admixture, we scan M\,M<i > 
100 GeV, \fj,\ < 1 TeV, and randomize the sign of fi, while keeping wia = mj = 3 TeV. The range of 
possible R are shown in Fig. gfor dark matter masses between 120 and 150 GeV. Points with wino 
fraction \Zy/\ 2 > 0.99 are plotted in black, points with Higgsino fraction \Zh u \ 2 + \Zn d \ 2 > 0.99 and 
tan/3 > 5 (tan/3 < 5) are plotted in blue (red), and points with max(<7 77 , l/2<r 7z ) > 10~ 32 cm 3 /s 
are plotted in gray. This lower bound on the cross section eliminates points that do not have a 
large enough cross section to explain the line; it is an extremely conservative choice, being roughly 
five orders of magnitude below the cross section obtained by |10j . 

Neutralinos that are dominantly wino or Higgsino give a sharp prediction for R th as a function 
of mass (for moderate to large tan/3). In particular, Higgsinos have a larger R t]l than winos because 
cr ann includes an 0(1) contribution from Z°Z° final states in addition to W + W~ , which is the only 
final state for the wino case. Note that the Higgsino-like prediction for R th is different than the 



one given in Eq. (12) due to the finite scan range for M^. Because the wino cross section to 7Z is 
larger than that for Higgsinos by two orders of magnitude, even a 1% wino admixture can change 
the annihilation cross section for the Higgsino-like points non-trivially. 

The large spread in R th for Higgsinos with small tan f3 arises from an interesting effect in the 
neutralino mass matrix. When tan/3 = 1, one of the Higgsino states is always pure and when 
\i < 0, this state is the lightest Higgsino. In this case, Mi can be small while still maintaining the 
purity of the Higgsino. Consequently, chargino and second neutralino mixing becomes non-trivial, 
lifting these masses above \x. The cross sections to 77 and jZ° are suppressed faster than those 
for annihilation to W + W~ and Z°Z°, which is why the red points are above the blue. Fig. [7] 
demonstrates that there is a lower bound on ii th in the decoupling limit with heavy sfermions, 
which is excluded by the constraint in Fig. |4j 

If we relax the assumptions on the masses of the sfermions and additional Higgs bosons, <r 77 
and a lz change at the percent level and <r ann is affected even less. Even when parameters are 
tuned so that the neutralino annihilation is dominated by the A resonance, the increase in the 
total annihilation cross section is orders of magnitude more then the increase in the rate for 77 or 
7Z . Hence, the lower bound on R th for dominantly wino/Higgsino neutralinos is robust. 

There is one final option to explore. When the sfermions are light, the bino is no longer inert. 
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FIG. 7: We have plotted R th 

wino fraction |Zn/| 2 > 0.99 are plotted in black, points with Higgsino fraction \Zh^, 
tan/3 > 5 (tan/3 < 5) are plotted in blue (red), and points with max(cr 77 , 1/2<t 7z 



Points with 



7z ) as a function of the neutralino mass. 

2 - 1 - 1 ZhA 2 > 0.99 and 



> 



10 -32 



cm' 



arc 



plotted in grey. We find a robust lower bound on cr ann /(2c7 77 + u lz ) as a function of mass for neutralinos 
with a large enough annihilation rate to 77 and/or jZ° to explain the 130 GeV line. For all points in this 
plot cr ann is dominated by some combination of a ww and cr zz ; the limits derived in Sec. Ill for \X ~^ W + W~ 
are relevant. 



In particular, for m x = 130 GeV and the slepton mass m~ f = 200 GeV, 7 



Bino: 



~ few x 10 -30 cm 3 /s; 
~ few x 10~ 31 cm 3 /s; 
~ a t fv ~ few x 10 -27 cm 3 /s. 



,R th ~ 10 3 . 



(13) 



The results of Figs. 10 and 11, which are relevant for lepton final states (see Appendix [C]) , exclude 
this possibility. However, the bino is not a good candidate for an even simpler reason. Recall that 
the value for cr 77 f found in [10J for the Einasto profile is 1.3 x 10 -27 cm 3 /s. Even allowing the 
large uncertainty in the shape of the profile at the Galactic Center, it is implausible that the cross 



sections given in Eq. (13) are large enough to yield the observation. As a result, the pure bino is 



not a good candidate to explain the 7-line. 

These arguments along with Fig. [7] demonstrate that there is a lower bound on a ww /(2 a. 
for neutralinos in the MSSM that have large enough cx 77 and/or <r 7z to be consistent with the data. 



77 a lz) 



Combined with the results of Sec. Ill, we see that the neutralino is excluded as an explanation of 
the 130 GeV Fermi line. 



V. CONCLUSIONS 

This paper presents constraints on the continuum photon spectrum for any dark matter can- 
didate used to explain the 130 GeV line in the Fermi data. Many models have been proposed 
to explain this feature |26[ l38H49j , but all must now satisfy the strong requirements on contin- 
uum annihilation derived here. In particular, we constrain the ratio of continuum photons versus 



Note we are neglecting the possibility of a Sommerfeld enhanced annihilation to photons, which can happen when 
the sfermion mass is very close to the neutralino mass. We thank Hai-Bo Yu for pointing out this possibility to us. 
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monochromatic line photons, a quantity that is independent of astrophysical uncertainties. The 
limit is strong enough to exclude models where the dark matter annihilates to W + W~ at tree level 
and into "yZ° and 77 at loop level. In particular, this rules out the neutralino explanation of the 
130 GeV 7 line. 

For W + W~ , Z°Z°, and bb annihilation, the data prefers no continuum contribution. For 
leptonic final states such as /^ + ^~ and t + t~, a small continuum contribution for a 145 GeV dark 
matter candidate provides the best fit to the data. Although the best fit models prefer little to no 
continuum, a model with line, background, and continuum contribution can still provide a good 
fit relative to a null power-law background. For example, the ratio of continuum photons from 
W + W~ to line photons can be as large as ~ 30 and still provide 4 a improvement over the null 
model for 130 or 145 GeV dark matter. A 3 a improvement over null is still achievable with a ratio 
about as large as 50. 

Figure [4] shows that the more aggressive shape constraint on continuum photons places a limit 
on the ratio of W + W~ , Z°Z° to 77 annihilations of O(10). Any model that is proposed to explain 
the 7 line must therefore have suppressed annihilation to these final states. Consider an example 
model where the dark matter has only one interaction — a Yukawa coupling to a new fermion / 
and scalar / that carry electroweak quantum numbers. If the mass of / and / are larger than the 
dark matter mass, the leading annihilation is loop suppressed. In particular, the 77, ^yZ°, Z° Z°, 
and W + W~ channels should all be 0(1). Although a more detailed computation is needed to 
determine the range of quantum numbers/couplings for / and /, this toy example illustrates some 
of the constraints that must be considered when model building for the Fermi line. 

There are reasons to be skeptical that this signal is the result of dark matter annihilations. The 
130 GeV 7 line was also observed along the Galactic plane along with evidence for other peaks 
in the gamma ray spectrum, implying that perhaps a more substantial look elsewhere correction 
should be applied [50]. Furthermore, a slight excess at 130 GeV was observed in data taken from 
the Earth's albedo, which could point to the presence of an unidentified instrumental bias for 
reconstructing photons at this energy [13]. The question remains whether the feature observed in 
|10j is a genuine dark matter signal, a systematic effect, or a new background |51j . 

If the 130 GeV line survives further scrutiny and is not a systematic or unknown background, 
then the statistical uncertainty on the line signal will be reduced as Fermi takes more data. This 
will tighten the constraint on the annihilation cross sections. In particular, the constraint on the 
ratio of 7Z to Z°Z° could approach 0(1). From the effective field theory point of view, it is 
difficult to understand how annihilations to photons will not also be accompanied by annihilations 
to Z° bosons. Hence, the updated versions of the constraints presented here could provide a 
significant challenge for model building a dark matter interpretation of the 130 GeV 7 line. 

Note Added 

At the time of completion of this manuscript, [52J appeared, which explores constraints on neu- 
tralino dark matter using the continuum spectrum from 20-200 GeV. Our results are in agreement. 
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Appendix A: Fermi Data Counts per Bin 
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TABLE II: Counts of Pass 7_Version 6 ULTRACLEAN events within 3° of the Galactic Center, binned in 
energy [GeV], with the inner 1° masked. Analysis cuts are as described in the text. 
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TABLE III: Counts of Pass 7_Version 6 ULTRACLEAN events within 3° of the Galactic Center, binned in 
energy [GeV]. Analysis cuts are as described in the text. 
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Appendix B: Modeling the Fermi Instrument Response Function 

FermVs Instrument Response Function describes the energy dispersion for a photon that hits 
the detector at a given energy and angle of incidence and can be modeled using the publicly 
available information in the Science Tools documentation. 8 The energy dispersion of the Fermi 
LAT is defined as 



5E _ E' — E 
J? ' E ' 



(Bl) 



where E is the true energy of the event and E' is the reconstructed energy. For the fits, a scaled 
dispersion x is used, where 



SE 



x 



E-S D (E,0) 



(B2) 



The scaling factor Sd depends on the true energy E and true incidence angle 6 of the event. It is 
fit with the function 



S D (E,e) = co\ log 



E 



[MeV] 



E E 

+ ci(cos 6) 2 + C2 log + C3 cos 6 + C4 log cos 9 + C5. (B3) 



[MeV] 



[MeV] 



For front-converted events, (c , c ± , c 2 , c 3 , c 4 , c 5 ) = (0.0210,0.0580,-0.207,-0.213,0.042,0.564). For 
back-converted events, (c , Ci, c 2 , c 3 , c 4 , c 5 ) = (0.0215, 0.0507, -0.220, -0.243, 0.065, 0.584). 

The energy dispersion depends on both the energy and incidence angle of the incoming photon. 
Its functional form is 



D ij (x) = { 



N l r(x,x 3 ,0£,7l) 
N t R (x,x J , of ,7/) 
N r R (x,x j , ov^Tr) 



if (x — xo) < —x 
if — x < (x — xo) < 
if < (x — xo) < x 
if (x — xq) > x 



(B4) 



where 



7 



R(x, xo, <t, 7) = — -, — ; - - — exp 
v u 2 1 /7r(l/ 7 ) C j ' 



X — Xo 

a 



(B5) 



The superscript i denotes the energy bin and j denotes the cos 6 bin. The values for x and the 
gammas are given by 



(x, 7l, li, 7r, 7ii) = (1-5, 0.6, 1.6, 1.6, 0.6). 



(B6) 



The normalization factors Nr., A;, N r , Nr are chosen such that D(x) is continuous and 
J_ D(x)dx = 1. The remaining five free parameters (xo, &l, <?h &r) are given by fits to the 
data for 18 separate energy bins and 8 separate cos# bins, for the front and back-converted events 
separately. 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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FIG. 8: The averaged energy dispersion D % (x) as a function of energy appropriate for a 130 GeV photon. 



The goal is to obtain the distribution for 5E, properly averaged over cos# as well as the front- 
and back-converted events. For a given energy bin, 



D 





Dl nt (x)X^ s> front + Dl ck (x)X^ SMck 




^ J ^efT,front + ^ J ^efT,back 





(B7) 



where N Cg is the number of bins in cos 9, X 3 is the exposure (summed over pixels within 3° of the 
Galactic Center), and A^ s is the effective area for the front and back converters of the tracker. 



Recall that x depends on cos 6 through Sd- Therefore, we substitute Eq. (B2 ) into the expression 
for D l (x), and evaluate So at the median value of cos# in the appropriate bin. If one specifies the 
true energy of the photon, E, where E is within the i th energy bin, then 



5E 



E ■ S D (E,cos9 med ) 
For reference, we plot the D l (x) for a 130 GeV photon in Fig. |8 



(B8) 
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Appendix C: Constraints for Alternate Final States 



XX bb 
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FIG. 9: The top plot gives the photon counts within 3° degrees of the Galactic Center with the inner degree 
masked. The solid red line shows the best fit model, which is given by the white cross in the bottom left 
plot. This best fit point has iV ann = 0; for reference the dashed black line shows the continuum spectrum for 
130 GeV dark matter annihilating into bb with an arbitrary normalization. On the bottom left, we show 1, 
2, and 3 a confidence regions (filled contours) for N ann / (N 17 + N lz ) as a function of mass for dark matter 
annihilation to bb. The ratio N 7z /N 17 is allowed to freely vary for each point in the grid. The solid black 
lines are the contours for N J7 +N J2 . . The best fit point is marked with a cross at m x = 130 GeV, # 7z / 77 = 0, 
and -/V ann = 0. On the bottom right, we show the shape analysis constraint. The shaded region corresponds 
to parameters where the fit is 2 a or worse with respect to the best fit point. 
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FIG. 10: The top plot gives the photon counts within 3° degrees of the Galactic Center with the inner 
degree masked. The solid red line shows the best fit model, which is given by the white cross in the bottom 
left plot. This best fit point has N ann ^ 0; the dashed black line shows the continuum spectrum for 145 
GeV dark matter annihilating into r + r _ with the best fit normalization. On the bottom left, we show 1, 
2, and 3 a confidence regions (filled contours) for N ann / (N 17 + N lz ) as a function of mass for dark matter 
annihilation to t + t~. The ratio N ryz /N 1J is allowed to freely vary for each point in the grid. The solid 
black lines are the contours for 7V 77 + N JZ . The best fit point is marked with a cross at m x = 145 GeV, 
6* 7z / 77 = 1.57, and N ann — 284. On the bottom right, we show the shape analysis constraint. The shaded 
region corresponds to parameters where the fit is 2 a or worse with respect to the best fit point. 
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FIG. 11: The top plot gives the photon counts within 3° degrees of the Galactic Center with the inner 
degree masked. The solid red line shows the best fit model, which is given by the white cross in the bottom 
left plot. This best fit point has iV ann ^ 0; the dashed black line shows the continuum spectrum for 145 
GeV dark matter annihilating into with the best fit normalization. On the bottom left, we show 1, 

2, and 3 a confidence regions (filled contours) for AT ann /(iV 77 + N lz ) as a function of mass for dark matter 
annihilation to The ratio N lz /N 11 is allowed to freely vary for each point in the grid. The solid 

black lines are the contours for N nl + N lz . The best fit point is marked with a cross at m x = 145 GeV, 
6* 7z / 77 = 1.57, and V ann = 160. On the bottom right, we show the shape analysis constraint. The shaded 
region corresponds to the parameters where the fit is 2 a or worse with respect to the best fit point. 
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